(N 

o 

(N 



Preprint typeset in JHBP style - HYPER VERSION OUTP-12-07P 



Gluon fusion contribution to W + W + jet production 



Tom Melia 

Rudolf Peierls Centre for Theoretical Physics, 1 Keble Road, University of Oxford, UK 



Email: b .melialOphysics .ox.ac.uk 



Kirill Melnikov 

Department of Physics and Astronomy, Johns Hopkins University, Baltimore, MD 
21218, USA 



Email: nelnikov@pha . jhu . edu 



^ \ Raoul Rontsch 

Rudolf Peierls Centre for Theoretical Physics, 1 Keble Road, University of Oxford, UK 



Email: |r . rontschl@physics . ox . ac . uk| 
Markus Schulze 

High Energy Physics Division, Argonne National Laboratory, Argonne, IL 604-39, USA 
| Email: [narkus . schulze@anl . gov 



Giulia Zanderighi 

' Rudolf Peierls Centre for Theoretical Physics, 1 Keble Road, University of Oxford, UK 

■ Email: |g . zanderighi l@physics . ox . ac .uk 

o 

Abstract: We describe the computation of the gg —> W + W g process that contributes 
to the production of two VF-bosons and a jet at the CERN Large Hadron Collider (LHC). 
While formally of next-to-next-to-leading order (NNLO) in QCD, this process can be eval- 
uated separately from the bulk of NNLO QCD corrections because it is finite and gauge- 
invariant. It is also enhanced by the large gluon flux and by selection cuts employed in 
the Higgs boson searches in the decay channel H — > W + W~, as was first pointed out by 
Binoth et al. in the context of gg —> W + W~ production. For cuts employed by the AT- 
LAS collaboration, we find that the gluon fusion contribution to pp — > W + W~j enhances 
the background by about ten percent and can lead to moderate distortions of kinematic 
distributions which are instrumental for the ongoing Higgs boson searches at the LHC. We 
also release a public code to compute the NLO QCD corrections to this process, in the 
form of an add-on to the package MCFM. 
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1. Introduction 

At the time of writing, the Large Hadron Collider (LHC) is running at a center-of-mass 
energy y/s = 8 TeV, having completed a very successful 7 TeV run at the end of 2011. 
Collected data are consistent with the Standard Model Higgs boson with a mass of 124—126 
GeV [1,2]. As more data are collected in 2012, the evidence for the new particle in this 
mass range will, hopefully, become stronger, opening up the way for a detailed exploration 
of its properties. 

Currently, the run = 125 GeV Higgs boson signal has the largest significance in 77 and 
ZZ* final states. However, eventually, the signal should also appear in the W + W~ final 
state thanks to the process gg —> H — )• W + W~ . In fact, understanding the properties of the 
new particle will require, at a minimum, the measurement of its branching fractions to as 
many final states as possible; therefore, measuring the branching fraction for H — > W + W~ 
is important. Higgs searches in the W + W~ mode are performed by binning final states 
according to the number of jets produced. This is primarily done for a detailed identification 
of backgrounds and for designing cuts that maximize the signal-to-background ratio for each 
of the jet bins. For zero- and one-jet bins, the dominant backgrounds are pp — > W + W~ 
and pp — > W + W~ + jet, respectively, while for the two-jet bin pp — > ti becomes dominant. 

QCD radiative corrections are known to change both signal and background rates by 
significant, jet-bin-dependent amounts. Higgs production rates are known through next- 
to-next-to-leading order (NNLO) QCD in the zero-jet bin [3-5] and through NLO QCD in 
the one- [6-8] and two-jet bins [9,10]. The production of a H^-pair without reconstructed 
jets is known through NLO QCD [11-15] and has been implemented in publicly available 
programs such as MCFM [16], MCONL0 [17], and POWHEG [18]. The production of a W-pair 
in association with one jet was calculated through NLO QCD in Refs. [19-21]. The NLO 
QCD corrections to the production of a VT-pair in association with two jets were computed 
in Refs. [22-24]. 
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Figure 1: Sample Feynman diagrams for gg — > W + W g through a fermion loop. 



Little is known about the main background processes beyond next-to-leading order. 
One exception is the study of the gluon-gluon fusion that gives rise to a pair of V7-bosons 
through a fermion loop [25-27]. Although this contribution appears at one-loop, it cannot 
contribute to next-to-leading order cross-sections because there is no tree process with two 
gluons in the initial state. Therefore, the gg — > W + W~ sub-process is finite and a gauge- 
invariant part of the NNLO contribution to W + W~~ production that can be computed 
independently. Thanks to the low production threshold of the two VF-bosons relative to the 
LHC center-of-mass energy, the cross-section arising from the gg — > W + W~ contribution 
is enhanced by the large gluon flux. To some extent, this may compensate the suppression 
by an additional power of a s and lead to the numerical importance of the gluon fusion 
sub-process. These NNLO corrections have been studied for the production of a W^-pair 
with no jets in Refs. [25-27]. It was found that the effect of these corrections is highly cut- 
dependent: for inclusive cuts they enhance the cross-section by 4 — 5%, which is comparable 
to the scale uncertainty at next-to-leading order, but the relative importance of the gluon 
fusion corrections increases dramatically if cuts designed to enhance the gg — >■ H —> W + W~ 
signal are applied. In fact, for the cuts of Refs. [26,27], gg — > W + W~ becomes the dominant 
radiative correction to pp —> W + W~ , overshooting the NLO QCD correction to this process 
by a factor of around seven. The gluon induced WW production is implemented in the 
generators MCFM [16] and gg2WW [27]; the latter can also be interfaced to parton shower 
programs. 

Since, as we pointed out earlier, the Higgs searches distinguish processes with differ- 
ent jet multiplicities and since the final state with a VF-pair and a jet is an important 
contributor to current Higgs searches, it is a relevant question if a similar enhancement of 
the gluon fusion contribution occurs there as well. The goal of this paper is to study this 
question. To this end, we compute the gluon-induced NNLO QCD corrections to 
jet production in proton collisions and explore its numerical significance in dependence of 
the applied cuts. 

To perform this calculation, we use a generalized unitarity implementation [28-31] of 
the Ossola-Pittau-Papadopoulos (OPP) procedure [32] for the reduction of tensor integrals 
(for a review see Ref. [31]). These on-shell techniques have had fantastic success over the 
past few years, leading to a large number of NLO QCD computations of very complex 
processes [33]. Recent developments in automating these and similar techniques for one- 
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A,(g v g 2 ,g 3 ) 



A,(g v g 3 ,g 2 ) 



Figure 2: Primitive amplitudes Af(gi, gi, g^) and Af(gi, 53, gz). One ordering of the VK-bosons 
is shown. However, we need to consider every insertion of the electroweak bosons relative to the 
ordered gluons. 



loop computations also seem very promising [34-37]. 

The remainder of this paper is organized as follows. In Section ^, we provide technical 
details of the calculation, pointing out a few subtleties that arise when dealing with a 
process that involves several colorless external bosons and virtual fermions. In Section ||, we 
present results for typical cuts that are used to observe weak gauge boson pair production 
and for a different set of cuts designed to increase the significance of the Higgs boson signal. 
We conclude in Section |j. In the Appendix, we provide numerical results for one-loop 
primitive amplitudes for further comparison. Finally, we note that as a by-product of this 
work we release a code that allows for the NLO QCD description of the pp — > W + W~+jet 
process . The code is in the form of a patch to the program MCFM [16], and can be 
obtained from [http:/ /www-thphys. physics.ox.ac.uk/people/TomMelia/tommelia.html , or 



from the MCFM website http:/ /mcfm. fnal.gov/ 



2. Technical Details 

We consider the gluon-gluon fusion contribution to the production of a PF-boson pair and a 
jet. At lowest order in QCD, this process occurs through a fermion loop and involves three 
gluons, see Fig. U This one-loop amplitude provides a finite, gauge-invariant contribution 
to the NNLO correction. We consider leptonic decays of W bosons and write the one-loop 
amplitude for the process — > gggW + {-^- v e e + )W~ (— > as 

1 



A){gi,g2,gz\v e ,e + ,n ,1^) = gl x 



(2.1) 

TT(T^T a ^)A f (g 1 ,g 2 ,g 3 ) + TT(T^T a ^)A f (g 1 ,g 3 ,g 2 )). 



Here g w and g s are the weak and strong couplings respectively, and the Af are color- 
ordered primitive amplitudes with all possible insertions of the two W bosons relative to 



1 Note that the public code does not include the NLO QCD corrections coming from fermion loops 
these contributions are several orders of magnitude smaller than the total cross-sections. 
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the ordered gluons, see Fig. |2|. The normalization of the SU(3) color group generators 
is Tr(T a T b ) = 5 ab . We take top and bottom quarks as massive, while the other four 
quarks are taken to be massless. We use a unit CKM matrix, so that each of the primitive 
amplitudes in Eq. ( |2.1[ ) contains two massless and one massive fermion loop, 

A f = 2A ffi + A fim . (2.2) 

It is also necessary to consider the production of WW through an intermediate Z-boson 
or a photon. Thus we can write the primitive amplitudes in Eq. (p.2[) as 



(2.3) 



A _ a IWW] f r [Z) . [Z V ] r [Z) . [Z A ] , rlA [7] 

A f,0- A f,ud + 2^ [ C V,q A f,q + C A,q A f,q + C q A f,q 
q£{u,d} 

A . - aWW] (c [Z] A [Zv] + C [z] A [Za] + C7*A M ' 

Af, m -A fitb + ^ \^V,q A f,q + ^A,q A f,q + U q A f,q 

q e{t,b} 

The first term refers to the amplitudes obtained by attaching the VF-bosons directly to 
the fermion loop. The remaining terms refer to amplitudes obtained for one flavor of 
quark circulating in the loop, with the PF-bosons attached through the vector and axial 
parts of a Z-boson or a photon, respectively. The couplings are Cy q = — 2Q g sin 2 9 W , 

C [Z l = -T g 3 , and Cp ] = 2Q q sin Qyji where Qq is the electromagnetic charge of the quark 
in the loop, 6 W is the weak mixing angle, and T^ ct = l/2,Tj sfe = —1/2. For the first 
term in Eq. fl2.3|) , it is necessary to consider every insertion of the VF-pair relative to the 
ordered gluons. There are twelve such insertions, allowing for the ordering W + W~ and 
W~W + . Similarly, for the remaining terms, we need to account for three orderings of 
an intermediate Z-boson or photon relative to the gluons. We account for the decays of 
W-bosons W(pi+pi) — > l(pi)l(pi) by constructing polarization vectors from lepton spinors. 
We write 

r(j>nPi) = „r" U, 2 \ .J?'"™ , (2-4) 



, n , _ -iu{pi)^{l-^)v{pj) 



where = {pi + pf) 2 is the invariant mass of the two leptons. The definitions of the 
polarization vectors of the intermediate Z-boson and the photon include a large number of 
terms that accommodate both single and double resonant amplitudes for the production 
of two VF-bosons. All amplitudes were checked against an independent home-made OPP- 
based program that can be used to compute Feynman diagrams directly. Furthermore, the 
term ^4j^j^ was checked against the publicly available program GoSam [34] for a few phase 
space points. 

We compute tree level amplitudes using Berends-Giele recursion relations [38]. There 
is a further technical aspect to be discussed, related to the use of these recursion relations 
with a cyclic amplitude. To illustrate this point, we consider two orderings of the W- 
bosons in the primitive amplitude Af(gi,g2,g3), with a different double cut on each, as 
shown in Fig. ||[ Both of these cuts split the one-loop amplitude into two tree-level helicity 
amplitudes, and we focus on the one involving g\ and W + . When a Berends-Giele current 
includes an electroweak boson, all insertions of that electroweak boson relative to the gluons 
are considered. In particular, this means that the tree-level amplitude arising from either 
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Figure 3: Two double cuts on a primitive amplitude with different orderings of the IT-bosons, 
leading to a potential source of double counting. 



cut shown will include the W + both before and after the gluon g%, reading clockwise. 
Since the amplitudes are cyclic, we would, in effect, be double counting. Because the 
cuts are on different propagators, it is difficult to identify automatically which cuts will 
end up duplicating one another in this manner. Instead, we observe that the cause of 
this problem is the cyclic nature of the amplitudes. However, looking at Eq. (2.1), it is 
clear that we can require gluon gi to be in a fixed position for both primitive amplitudes, 
with appropriate ordering of the other two gluons and all insertions of the VT-bosons. We 
modify our Berends-Giele currents accordingly, so that the positions of the electroweak 
bosons are fixed relative to g±, but still allowing all insertions relative to the other two 
gluons. Referring again to Fig. [3L the left amplitude will have the Ty + -boson only before 
the gluon (again, reading clockwise), while the right amplitude will have the T^ + -boson 
only after the gluon, and the double counting is avoided. 



3. Results 



In this Section we present the results for the process pp — > W + W~ + jet at the LHC, 
including the gluon fusion contribution gg — > W + W~g. We consider two center-of-mass 
energies, -y/s = 8 GeV and ^/s = 14 TeV and let the VF-bosons decay into a pair of 
leptons of definite flavor, W + W~ — > u e e + pTv^. We note that the cross-section inclusive in 
lepton flavors = {e + e~ ,e + fi~ , u, + e~ , n, + fi~} is a factor of four larger than the results 
presented below. We take the VF-boson mass and width to be mw = 80.399 GeV and 
T\y = 2.085 GeV, while the Z-boson mass and width are taken to be mz = 91.1876 GeV 
and Tz = 2.4952 GeV. We use a top quark mass mt = 172.9 GeV and a bottom quark 
mass nib = 4.19 GeV and take all other quarks to be massless. The weak couplings are 
defined through the Fermi constant G F = 1.166364 x 10" 5 GeV~ 2 with gl = %G F m^/\/2, 
and the weak mixing angle is given by sin 2 8 W = 1 — m^/m^. Jets are defined through 
the anti-/c^ algorithm [39] with R = 0.4, as implemented in Fast Jet [40]. We use MSTW08 
parton distribution set [41] and we employ leading, next-to-leading and next-to-next-to- 
leading order parton distribution functions (pdf) to compute the relevant contributions to 
the cross-sections. We stress that, for the purposes of this discussion, the gluon fusion 
process gg — > W + W~g is a NNLO contribution and therefore it is computed with NNLO 
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Standard Cuts 









incl 
a NLO 


cxcl 
a NLO 


<5cnnlo 


^NNLO/^NLO 


8 TeV 


WW 
WW] 


141.0(1)+!* 
87.8(1);™ 


232.0(4);^ 
111.3(2)"^ 


143.8(2)+^ 
66.6(2)1^ 


3-4(l)+l;° 


3.5% 
3.1% 


14 TeV 


WW 
WW] 


259.6(2)™ 
203.4(1)7™ 


448.3(5)7™ 6 
254.5(4);i° 2 


242.0(3)1^ 
127.6(4)™ 


23-6(1);^ 

n.8(4);f 7 


5.3% 
4.6% 



Table 1: Cross-sections for pp — > W + (v e e + )W~(fj,~i>fj,) + n jets, n = 0, 1 at the y/s = 8 TeV 
and yfs = 14 TeV LHC, using the standard cuts described in the text. All cross-section values 
are in femtobarns. The central values are computed with factorization and renormalization scales 
fi = 2m\v, with the statistical error shown in parentheses. The upper and lower values for cross- 
sections are obtained with the choice of scale fi = my/ (subscript) and /i = Amy/ (superscript). The 
last column shows the relative size of the NNLO contribution and the inclusive NLO cross-section 
for the central scale choice. 

pdfs 2 . 

We first consider cuts that are typical for processes with W-bosons and jets, and refer 
to them as "standard cuts". We require that i) both leptons have transverse momenta 
Pxi > 20 GeV and pseudorapidity \r/i\ < 2.5, ii) the missing transverse momentum satisfies 
PT,miss > 30 GeV and Hi) jets have transverse momenta pr,j > 20 GeV and pseudorapidity 
\r)j\ < 3.2. In Table |l|, we show the LO and NLO cross-sections for pp — > W + W~+ jet, as 
well as the gluon fusion (NNLO) contribution, using these cuts for two LHC energies. Also 
shown are the sizes of the NNLO corrections relative to the inclusive NLO cross-sections. 
We set the factorization and renormalization scales equal to one another, [ip = /i/j = /x, 
and show the scale variation miy < fj, < Amy/ with a central scale choice n = 2miy. We 
also show the cross-sections for pp — > W + W~ without jets in Table |l[ for comparison 
with the earlier calculation of Ref. [26]. Although the setup here is slightly different, in 
particular cuts on missing transverse energy, parton distribution functions and some input 
parameters, our results are similar to the ones reported in Ref. [26] for LO, NLO and the 
gluon fusion contributions to pp —> W + W~ process. It follows from Table |j] that the gluon 
fusion corrections are modest; they change the NLO cross-sections by 3 — 5% which is to be 
compared with 0(65 — 70%) inclusive NLO QCD corrections in the case of pp — > W + W~ 
and 0(15 — 25%) in the case of pp — > W + W~ + jet. The relative importance of the 
gluon fusion contribution grows with the energy which is to be expected because of the 
rapid growth of the gluon flux at small values of x ~ m w + w - / yfs. For both processes, 
pp — > W + W~ and pp — > W + W~ + j, the gluon fusion corrections are comparable to 
the scale uncertainty of the NLO cross-section at ^/s = 8 TeV and slightly exceed it at 
sfs = 14 TeV. 

Next, we consider cuts similar to those used by the ATLAS collaboration in their 

2 We note that the use of NNLO pdfs for the computation of the gluon fusion contribution leads to a 
smaller result (by roughly 30%) compared to that obtained using LO pdfs, since the NNLO gluon flux is 
smaller than the LO flux. We nevertheless use NNLO pdfs since this is what would be done once the NNLO 
description of the W + W~ production becomes available. 
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Higgs search in the W + W channel [42] ("Higgs search cuts"). To this end, we require 
that i) the hardest lepton has Pt,i, max > 25 GeV, the softest lepton has pr,l,mia 

> 15 GeV, 

ii) the lepton rapidities are \rji\ < 2.5 3 ; Hi) the relative azimuthal angle between the 
two leptons is small A(ftu < 1.8, iv) the two leptons are separated from each other by 
Ai?;z = yj A^| + A 77^ > 0.3, v) the invariant mass of the charged lepton system must satisfy 
10 GeV < mu < 50 GeV, vi) the missing relative transverse energy satisfies £K > 25 
GeV. We note that the missing relative transverse energy is defined as 

^l = |p miss |sinA«/» min , (3.1) 

where A0 mm = min(A(/) m ; ss , 7r/2) and A0 m i ss is the azimuthal angle between the missing 
transverse momentum vector p™ lss and the momentum of the nearest charged lepton or 
jet with pt > 25 GeV. Jets are identified using the anti-/cj_ algorithm with R = 0.4 and 
Ptj > 25 GeV. All jets must have pseudorapidities in the region \rjj\ < 4.5. For the WW 
process we require additionally that the charged lepton system has a transverse momentum 
p T ,ii > 30 GeV. 

The cross-sections obtained using these cuts are shown in Table ^. While the impor- 
tance of NLO QCD corrections is considerably reduced with these cuts, the relative impor- 
tance of gluon fusion corrections increases by a factor of two, compared to the standard cuts. 
The change to the exclusive NLO cross-sections for pp — > W + W~ and pp — > W + W~ + j 
because of the gluon fusion contribution is about 6% for \fs = 8 TeV and is close to 10% 
for y/s = 14 TeV. Therefore, for the Higgs search cuts, gluon fusion corrections are larger 
than the scale variations of the exclusive NLO QCD cross-sections for the process 
and comparable to the scale variations of the exclusive NLO QCD cross-sections for the 
process W + W~ + 1 jet. 

Furthermore, using the Higgs search cuts described above and a Higgs mass of mu = 
125 GeV, the signal H — > WW cross-section (in the e + fi~ channel considered here) at 
NLO QCD is approximately 5 fb (12 fb) at 8 TeV (14 TeV) in the 0-jet channel and 2 fb 
(5 fb) at 8 TeV (14 TeV) in the 1-jet channel. The gluon fusion corrections to the WW 
background alone therefore amount to almost half of the signal cross-sections. 

We note that the gluon fusion corrections to pp — > W + W~ that we find in this paper 
are significantly smaller than O(30%) effects reported in Refs. [26, 27]. This happens 
because the cuts used in those references are more aggressive than the cuts that we use 
here. In particular, stronger cuts on A0//, mu and, especially, on the transverse momenta 
of the charged leptons, account for the larger gg — > WW contribution. We have checked 
that with cuts similar to those adopted in Ref. [26, 27] we indeed would have a significant 
increase in the relative importance of the gluon fusion contribution. 

Finally, the integrated cross-section - even when computed with cuts optimized for 
the Higgs boson searches - is not the whole story since experimentalists rely on shapes 

In the measurements of ATLAS [42] the rapidity region 1.37 < jr;| < 1.52 is excluded for electrons. For 
simplicity we use instead full acceptance both for electrons and muons. We also ignore issues related to 
specifics of the detector. Our treatment of the lepton separation is also slightly simpler than the one in 
Ref. [42]. 
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Figure 4: Distributions of the mass of the charged lepton system mu, transverse mass of the 
T^-pair m,T,ww, the azimuthal angle between the leptons Ac/);/ and the transverse momentum of 
the positron p T e + , shown at NLO accuracy with and without the fermion loop NNLO contribution. 
We use the Higgs search cuts as described in the text, and we display results for the 8 TeV LHC. 
The upper and lower bands show the maximum and minimum deviations from the central scale 
value 2mw- The ratio -&Tnlo as defined in Eq.(3.3) is also displayed. 



Higgs search cuts 







0"LO 


incl 
NLO 


excl 
"NLO 


<5tTNNLO 


5cnnlo/o"nlo 


8 TeV 


WW 
WW] 


35.6(1)1?! 
12.6(1)-};! 


51.1(1)^:9 

10.8(1)1$ 


38.8(1)1^ 
10-6(1)1^ 


2.7(1)^; 5 7 

0-6(1);°,;! 


7.0% 
5.7% 


14 TeV 


WW 
WWj 


63.4(1)11;? 
28.7(1)+1;| 


91.9(2);°i 
21.6(1)1*;? 


63.4(2)±1;5 
20.5(1)1 1 ! ; 2 7 


7.5(1);^ 

1.8(2)^7 


11.8% 
8.8% 



Table 2: As for table 1, but using the ATLAS Higgs search cuts as described in the text. The 
ratio in the last column is taken using the exclusive NLO cross-section. 



of distributions to distinguish signal and backgrounds. In Fig. || we show the relevant 
kinematic distributions for the background W + W~ +j process, paying particular attention 
to the modifications caused by the gluon fusion contribution. We display distributions in 
the azimuthal angle between the charged leptons, the invariant mass of the charged leptons, 
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the transverse momentum of the positron, and the transverse mass of the W + W pair 
defined as 



mr,WW = \ (Et,U + PT.miss) 2 - |PT,miss + Pt,z/| 2 - (3.2) 



Here E^u = y\PT,u\ 2 + m u an dpr,miss = |pT,mlss|- Since spin correlations in H — > W + W~ 
decays force the two charged leptons to be close to each other, the most interesting, signal- 
rich regions of the A(f>u- and m^-distributions are at small values of the corresponding 
variables. The distribution in the transverse mass of the W + W~ pair is also useful be- 
cause restricting this variable to values comparable to the Higgs boson mass removes the 
interference between the Higgs signal and the background from direct W + W~ produc- 
tion [43]. Furthermore, ATLAS collaboration uses the distribution of rriT,ww to test for 
the presence of a signal [42] . In the lower panes in Fig. ^j, we show the differential AT-factors 
defined as 

do"NLO+(5NNLO 



K 



NLO 



d<7NLO 



(3.3) 

fi=2mw 



We notice that the ATnlo factor spread over the distributions mu, A(f>u and px e + 1S close to 
uniform, while the relative importance of the 5o"nnlo contribution increases with m,T t ww- 



4. Conclusion 

In this paper, we considered the gluon fusion contribution to the production of a pair of 
VF-bosons and a jet at the LHC. While formally part of the NNLO QCD correction to 
the pp — > W + W~+jet process, this contribution can be treated separately because it is 
finite and gauge invariant. In addition, gluon fusion contributions are enhanced by the 
large gluon flux at the LHC and may therefore be an important part of the backgrounds to 
Higgs boson searches. In fact, it was found in Refs. [26,27] that, for a particular set of cuts, 
the gg — > W + W~ contribution becomes the largest radiative correction to pp — > W + W~ , 
changing the leading order contribution by as much as 30%. 

In the current paper, we used the set of cuts employed by the ATLAS collaboration 
in their current searches for the Higgs boson to estimate the gluon fusion contribution 
to pp —> W + W~ + or 1 jet and we find significantly smaller radiative corrections than 
the results quoted in Refs. [26,27]. This is a combined effect of looser cuts on Acftu, mu 
and the lepton transverse momenta currently used by the ATLAS collaboration compared 
to those used for the discussion in Refs. [26,27]. For the ATLAS Higgs search cuts we 
find that the gluon fusion contribution to pp — > W + W~j is somewhat larger than, but 
still comparable to, the scale uncertainty of the exclusive NLO cross-section and most of 
the kinematic distributions. However, as our results show, these conclusions very much 
depend on the exact cuts applied. Should the experimental cuts change significantly, the 
calculations reported in this paper should be reconsidered. 
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the importance of gluon fusion background processes for the Higgs boson searches. This 
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06CD11357, by the LHCPhenoNet network under the Grant Agreement PITN-GA-2010- 
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\Af,o\ 


1 Af,m 1 


9i,92,93 
9i,92i93 
9i,92~,93 


3.844351 
5.699256 
3.522647 
23.21330 


0.9542655 
1.606166 
6.204687 
26.65234 



Table 3: Numerical results for the massless and massive fermion loop primitive amplitudes 
Af,o{gi,g2,g3) and Af, m (gi,g2,g3,), in units of 10~ 6 GeV~ 3 . 



264564, and by the British Science and Technology Facilities Council. M.S. is grateful for 
support from the Director's Fellowship of Argonne National Laboratory. 

A. Appendix 

In this appendix, we display the results for some of the primitive amplitudes, as well as 
a full amplitude summed over helicities, for a single phase space point. We consider the 
process gg — > W + (— > u e e + )W~ (— > fJ-~v^)g, and use the phase-space point (E,p x ,p y ,p z ) 
(in GeV): 

p Bl = (-500,0,0,-500), 
Pg2 = (-500,0,0,500), 

Pg3 = (86.3540681437814, -15.2133893202618,37.6335512949163, -76.2187226821854), 
p Ve = (414.130068374543, 232.145564945939, 332.7544367808, -82.9857518524426), (A.l) 
Pe+ = (91.8751521026384, -43.3570226791011, -24.0058236140057, 77.3623460793435), 
Pp- = (280.118181809376, -83.1261116505822, -263.203856758651, 47.7490851160266), 
p Pfi = (127.522529569661, -90.4490412959935, -83.1783077030789, 34.0930433392580). 

We present the results for the primitive amplitude Af(gi,g2,gs) as defined in Eq. 
for four sets of helicities of the three gluons. In Table ||, we display both massless and 
massive fermion loop amplitudes as defined in Eq. fl2.2j). We show only the finite parts of 
the amplitude, since there are no poles. 

We also give results for the full amplitude squared, that is, the color-squared amplitude 
from Eq. (^4]) summed over all gluon helicities. If we consider only one massless generation 
in the loop, the result is 8.348897 x 10~ 8 GeV~ 6 . If we consider one massive generation 
in the loop, the result is 1.128414 x 10~ 7 GeV~ 6 . For the setup used in deriving our 
results, i.e. two massless and one massive generation in the loop, the amplitude squared is 
7.830968 x 10- 7 GeV~ 6 . 
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